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Abstract. The closed time-path formalism is a powerful Green's function formulation 
to describe non-equilibrium phenomena in field theory and it leads to a complete non- 
equilibrium quantum kinetic theory. We make use of this formalism to write down the 
set of quantum Boltzmann equations relevant for leptogenesis. They manifest memory 
effects and off-shell corrections. In particular, memory effects lead to a time-dependent 
CP asymmetry whose value at a given instant of time depends upon the previous history 
of the system. This result is particularly relevant when the asymmetry is generated by 
the decays of nearly mass-degenerate heavy states, as in resonant or soft leptogenesis. 



1. Introduction 

Thermal leptogenesis [H El El H] is a well-motivated scenario for the production of the 
baryon asymmetry in the early Universe. It takes place through the decay of heavy 
right-handed (RH) Majorana neutrinos. The out-of-equilibrium decays occur violating 
lepton number and CP, thus satisfying Sakharov's conditions [5]. In grand unified theories 
(GUT) the masses of the heavy Majorana neutrinos masses are typically smaller than the 
scale of unification of the electroweak and strong interactions, Mqut ~ 10 16 GeV, by a 
few to several orders of magnitude. This range coincides with the range of values of the 
heavy Majorana neutrino masses required for a successful thermal leptogenesis. Being RH 
neutrinos a key ingredient in the formulation of the well-known seesaw mechanism [6], 
which explains why neutrinos are massive and mix among each other and why they 
turn out to be much lighter than the other known fermions of the Standard Model 
(SM), thermal leptogenesis has been the subject of intense research activity in the last 
few years. For instance, flavour effects have been recently investigated in detail (and 
shown to be relevant) [Zl[Sl[Sl[ini[IIl[I2l[I3l[Ill[ia[IQ[IIl[ia[iniffl including the 
quantum oscillations/correlations of the asymmetries in lepton flavour space 0,[IE]- The 
interactions related to the charged Yukawa couplings enter in the dynamics by inducing 
nonvanishing quantum oscillations among the lepton asymmetries in flavour space [9]. 
Therefore the lepton asymmetries must be represented as a matrix in flavour space, the 
diagonal elements are the flavour asymmetries, and the off-diagonals encode the quantum 
correlations. The off-diagonals should decay away when the charged Yukawa couplings 
start mediating very fast processes. For instance, at temperatures T ~ 10 12 GeV, the 
interactions mediated by the tau- Yukawa coupling enter in equilibrium and the tau-flavour 
becomes distinguishable from the muon and the electron flavour. A full treatment of this 
transition based on the quantum, rather than classical, Boltzmann equations is suitable 
to properly describe all the physical effects. 

Mainly, but not only, motivated by the impact of flavour effects, in this paper we 
set the stage for the study of the dynamics of thermal leptogenesis by means of quantum 
Boltzmann equations (for a previous study, see [ZlJ). They are obtained starting from the 
non-equilibrium quantum field theory based on the Closed Time-Path (CTP) formulation. 
We will see that the resulting kinetic equations describing the evolution of the lepton 
asymmetry and the RH neutrinos are non-Markovian and present memory effects. In 
other words, differently from the classical approach where every scattering in the plasma 
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is independent from the previous one, the particle abundances at a given time depend 
upon the history of the system. The more familiar energy-conserving delta functions are 
replaced by retarded time integrals of time-dependent kernels and cosine functions whose 
arguments are the energy involved in the various processes. Therefore, the non-Markovian 
kinetic equations include the contribution of coherent processes throughout the history of 
the kernels. 

We will see that one immediate consequence of the CTP approach to thermal 
leptogenesis is that the CP asymmetry is a function of time and its value at a given 
instant depends upon the previous history of the system. Furthermore, in the quantum 
approach, the relaxation times are expected to be typically longer than the one dictated 
by the classical approach. Particle number densities are replaced by Green's functions 
which are subject both to exponential decays and to an oscillatory behaviour. This 
restricts the range of time integration for the scattering terms and leads to a decrease of 
the wash-out rates. This is a well-established fact in nuclear collisions |22j. If the time 
range of the kernels are shorter than the relaxation time of the particles abundances, 
the solutions to the quantum and the classical Boltzmann equations differ only by terms 
of the order of the ratio of the timescale of the kernel to the relaxation timescale of 
the distribution. In thermal leptogenesis this is typically the case. However, there are 
situations where this does not happen. For instance, in the case of resonant leptogenesis 
and soft leptogenesis, two RH (s)neutrinos Ni and N 2 are almost degenerate in mass and 
the CP asymmetry from the decay of the first RH Ni is resonantly enhanced if the mass 
difference AM = (M 2 — Mi) is of the order of the decay rate of the second RH neutrino 
rjv 2 . The typical timescale to build up coherently the CP asymmetry is of the order of 
1/AM, which can be larger than the timescale ~ l/IVi for the change of the abundance 
of the Nis. This tells us that the reduction of the quantum Boltzmann equations to the 
classical ones should not be taken for granted. 

The paper is organized as follows. In Section [2] we provide the basics of the CTP 
formulation of non-equilibrium quantum field theory, we set the notations and state the 
main results which will be used throughout the paper. The quantum Boltzmann equation 
describing the time evolution of the lepton number asymmetry is formally derived in 
Section [3j in terms of propagators and self-energies of the leptons. Sections H] and [5] 
are devoted to write down more explicitly the quantum Boltzmann equations for the 
abundance of right-handed neutrinos and the lepton asymmetry, respectively. The familiar 
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classical equations are then easily recovered as a limit case. Section also contains 
an analysis of the time- dependent CP asymmetry obtained in this more general setup. 
Concluding remarks are in Section [6j 



2. The CTP formalism for non-equilibrium quantum field theory 



We first briefly present some of the basic features of the non-equilibrium quantum field 
theory based on the CTP formalism, also known as Schwinger-Keldysh formalism [23] . 
The interested reader is referred to the excellent review by Chou et al. [21] for a more 
exhaustive discussion. Since we need the time evolution of the particle asymmetries with 
definite initial conditions and not simply the transition amplitude of particle reactions, the 
ordinary equilibrium quantum field theory at finite temperature is not the appropriate 
tool. The most appropriate extension of the field theory to deal with non-equilibrium 
phenomena amounts to generalizing the time contour of integration to a closed time- 
path. More precisely, the time integration contour is deformed to run from to to +oo and 
back to to (see Fig. [1]). The CTP formalism is a powerful Green's function formulation 
for describing non-equilibrium phenomena in field theory. It allows to describe phase- 
transition phenomena and to obtain a self-consistent set of quantum Boltzmann equations. 
The formalism yields various quantum averages of operators evaluated in the in-state 
without specifying the out-state. On the contrary, the ordinary quantum field theory 
yields quantum averages of the operators evaluated with an in-state at one end and an 
out-state at the other. 

For example, because of the time-contour deformation, the partition function in the 
in-in formalism for a complex scalar field is defined to be 



Z[J,J*] =Tr 
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where C in the integral denotes that the time integration contour runs from to to plus 
infinity and then back to to again. The symbol p represents the initial density matrix and 
the fields are in the Heisenberg picture and defined on this closed time-contour (plus and 
minus subscripts refer to the positive and negative directional branches of the time path, 
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Figure 1. The complex time contour for the CTP formalism. 

respectively). The time-ordering operator along the path is the standard one (1+) on the 
positive branch, and the anti-time-ordering (71) on the negative branch. As with the 
Euclidean-time formulation, scalar (fermionic) fields are still periodic (anti-periodic) 
in time, but with 4>(t,x) = <p(t — i/3,x), [3 = 1/T. The temperature T appears due to 
boundary condition, and time is now explicitly present in the integration contour. 

We must now identify field variables with arguments on the positive or negative 
directional branches of the time path. This doubling of field variables leads to six different 
real-time propagators on the contour [24J. These six propagators are not independent, 
but using all of them simplifies the notation. For a generic charged scalar field they are 
defined as 

G>(x,y) = -Gj + {x,y) = -i(<l>(x)<jJ{y)), 

G<(x,y) = -G+-{x,y) = -i(tf(y)<l>(x)), 

Gl(x, y) = G++ (x, y) = 6(x, y)G^(x, y) + 9(y, x)Gf{x, y), 

Gl(x, y) = G~f {x, y) = 6{y, x)G>{x, y) + 6{x, y)G<(x, y), 

G r ^{x, y) = G^ — G^ = G^ — G^, G^(x, y) = G^ — G^ = G^ — G^, (2) 

where the last two Green's functions are the retarded and advanced Green's functions 
respectively and 8(x, y) = 9(t x — t y ) is the step function. 

For a generic fermion field ip the six different propagators are analogously defined as 

G ^(x,y) = -G^ + (x,y) = -i($(x)$(y)) , 

G i>( x ,y) = - G y (x,y) = +i(ip(y)ilj(x)}, 

G%{x, y) = Gj + (x, y) = 6(x, y)G^{x, y) + 0(y, x)G^(x, y), 

G%{x, y) = Gy (x, y) = 0(y, x)G^(x, y) + 9(x, y)G$(x, y), 

G^,(x, y) = G^ — G^ = G^ — G^, G^(x, y) = — G^ = G^ — G^. (3) 

From the definitions of the Green's functions, one can see that the hermiticity properties 

(ij G,p(x,y)) = ij°G^(y,x), (iG^(x,y))' f = iG^(y,x). (4) 
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are satisfied. When computing a loop diagram, one has to assign to the interaction points 
a plus or a minus sign in all possible manners and sum all the possible diagrams, taking 
into account that, by definition, vertices with a minus sign must be multiplied by —1. 

For equilibrium phenomena, the brackets (• • •) imply a thermodynamic average over 
all the possible states of the system. For homogeneous systems in equilibrium, the Green's 
functions depend only upon the difference of their arguments (x, y) = (x — y) and there is 
no dependence upon (x+y); for systems out of equilibrium, the definitions and ([3]) have 
a different meaning. The concept of thermodynamic averaging is now ill-defined. Instead, 
the brackets mean the need to average over all the available states of the system for the 
non-equilibrium distributions. Furthermore, the arguments of the Green's functions (x, y) 
are not usually given as the difference (x — y). For example, non-equilibrium could be 
caused by transients which make the Green's functions depend upon (t x ,t y ) rather than 

(t x — ty). 

For interacting systems, whether in equilibrium or not, one must define and calculate 
self-energy functions. Again, there are six of them: £*, £*, S < , S > , Y7 and E a . The same 
relationships exist among them as for the Green's functions in (|2j) and ([3]), such as 

S r = E* - £< = £> - £*, S a = S* - £> = £< - £*. (5) 

The self-energies are incorporated into the Green's functions through the use of Dyson's 
equations. A useful notation may be introduced which expresses four of the six Green's 
functions as the elements of two-by-two matrices 




G= \ £ % ), E=|;s ~* h (6) 

where the upper signs refer to the bosonic case and the lower signs to the fermionic case. 
For systems either in equilibrium or in non-equilibrium, Dyson's equation is most easily 
expressed by using the matrix notation 

G(x,y) = G°(x,y) + [ d 4 z, [ d 4 z 2 G°(x, Zl )%( Zl , z 2 )G(z 2 ,y), (7) 



where the superscript "0" on the Green's functions means to use those for noninteracting 
system. It is useful to notice that Dyson's equation can be written in an alternative form, 
instead of ([7j), with G° on the right in the interaction terms, 

G(x,y) = G°(x,y) + J d 4 z 3 J d 4 z 4 G(x,z 3 )Z(z 3 ,z,)G\z 4 ,y). (8) 



Eqs. ((?]) and (|S]) are the starting points to derive the quantum Boltzmann equations 
describing the time evolution of the CP-violating particle density asymmetries. 

For a generic complex scalar field we will adopt the real-time propagator in the form 
Crl(k, t x — t y ) in terms of the spectral function p^(k, fc ) [21] (for the equilibrium analog, 
see [26J) 

/+°° ,71,0 
-oo 2n 

X { [1 + U(k )} 6{t x - ty) + U(k )6(ty ~t X )}, (9) 

where f<j,(k°) represents the scalar distribution function. Notice that the expression Q) 
is only strictly correct for free Green's functions. To account for interactions with the 
surrounding particles of the thermal bath, particles must be substituted by quasiparticles, 
dressed propagators are to be adopted and self-energy corrections to the propagator 
modify the dispersion relations by introducing a finite width T^k). In the limit of small 
decay width, the spectral function is expressed by 



where u%(k) = k 2 + M|(T) and Afy(T) is the thermal mass. Performing the integration 
over k° and picking up the poles of the spectral function (which is valid for quasi-particles 
in equilibrium or very close to equilibrium), one gets 

G>(k,t x -t y ) = {[1 + Ufa - iT+)] e -**-«W-«»> 

+ U{u lf> + iT 4> )e i ^ +ir ^-^}, 
G<(k,t y -t x ) = -±-{Ufo + t Tt) e -^m^-ty) 

+ [1 + Ufa-iT^)] e *K-NT*)(*--*»)}, (11) 

where and denote the distribution function of the scalar particles and antiparticles, 
respectively. The expressions ( TTTT) are valid for t x — t y > 0. 

Similarly, for a generic fermion ip, we adopt the real-time propagator in the form 
Gv,(k, t x — t y ) in terms of the spectral function p^(k, k ) [21] (for the equilibrium analog, 
see [26J) 

^-e-^'^p^k ) 

X { [1 - U(k )] 6(t X - ty) - U(k )6(ty ~t X )}, (12) 
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p (k, k°) = i 



where f^{k°) represents the fermion distribution function. Again, particles must be 
substituted by quasiparticles, dressed propagators are to be adopted and self-energy 
corrections to the propagator modify the dispersion relations by introducing a finite width 
Tjp(k). Equilibrium fermionic dispersion relations are highly nontrivial [27] and here we 
will adopt relatively simple expressions for the fermionic spectral functions holding in the 
limit in which the damping rate is smaller than the real part of the self-energy of the 
fermionic excitation [28]. For a fermion with chiral mass m^, it reads 



p v ,(k, k°) = i($ + 7710) 



1 1 



(fcO + %e + lT ^2 _ u 2^ k) (£0 _ le _ lT ^2 _ 



, (13) 



where uj^(k) = k 2 + Mj(T) and M^(T) is the effective thermal mass of the fermion in the 
plasma (not a chiral mass). We reiterate the fact that at finite temperature the dispersion 
relation is in fact more complicated than f[T5j) . However, the latter suffices for our goals. 
Performing the integration over k° and picking up the poles of the spectral function (which 
is valid for quasi-particles in equilibrium or very close to equilibrium), one gets 

G>(k,t x -t y ) = -J-{^ + m ^)[l-f^-tT^] e -K-«»(^) 

+ 7° - 77^) 7 %K + e^ +lT ^-**)} , 
G<(k, ty -t x ) = J-{(t + n-ty) + il» e-^-W-M 

+ 7° (# - m*) 7° [1 - ~ e^*^)^-**)} , (14) 

where k° = u>^ and fy, denote the distribution function of the fermion particles and 
antiparticles, respectively. The expressions (fl4"|) are valid for t x — t y > 0. 

The above definitions hold for the Higgs and lepton doublets (after inserting the 
chiral left-handed projector Pi), as well as for the Majorana RH neutrinos, for which 
one has to assume identical particle and antiparticle distribution functions and insert the 
inverse of the charge conjugation matrix C in the dispersion relation. 



3. The quantum Boltzmann equation for the lepton asymmetry: generalities 

Since we are interested in the lepton asymmetry, we define the fermionic CP-violating 
current of the lepton doublet i{ as 



(j£{z)) = (£ t (x)Y^(x)} = (n Ci (x),J Ci (x))- (15) 
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The zero-component of this current represents the number density of particles minus 
the number density of antiparticles and is therefore the relevant quantity for leptogenesis. 

We want to find a couple of equations of motion for the interacting fermionic Green's 
function when the system is not in equilibrium. Such equations may be found 

by applying the operators i $ x and i $ y on both sides of Eqs. ((Tj) and (jSJ), respectively. 
We find 

i$ x G £i (x,y) = 5 {4) (x,y)I + J c^zE^x, z)G ei (z,y), 

Gt.{x,y)i $ y = - 5 {4 \x,y)I - J d 4 zG £i (x, z)E ei {z, y), (16) 

where I is the 4x4 identity matrix. We can now take the trace over the spinorial indices 
of both sides of the equations, sum up the two equations above and finally extract the 
equation of motion for the Green's function Gf. 



Tr 



i $ x +i 



Gl(x,y) 



(TzTr 



E>(x,z)G t £i (z,y)-El i (x,z)Gl(z,y) 



(17) 



- G>(x,z)Zl(z,y) + G%(x,z)Z>(zy : y) . 

Since we want to compute the average of observables at equal time, we will identify the 
variables x and y. Therefore, it turns out useful to define a center-of-mass coordinate 
system 

X = (t,X) = ^(x + y), (T,f) = x-y. (18) 
The notation of the Green's function is altered in these center-of-mass coordinates 



G/U.u) G, [r.F.t.X] [l- ] - T ,X- ] -f 



t+ l -r,X+ l -r)). (If)) 



Making use of the center-of-mass coordinate system, we can work out the left-hand side 
of Eq. (TTJ 



Tr 



G>(t,X,r,r) + G>(t,X,r,r)i $ y 

d 



i{dF + dl)i{i il n i )\ T __ 



T=f=0 

(itW-fiiix)) - 



•Jo. ■ 



(20) 



=f=o qx»" y " y " dX» 
The next step is to employ the definitions in ([3]) to express the time-ordered functions 
G\., G\., S^., and H\. in terms of Gf., Gf., Tlf. and Gf.. The time integrals are separated 
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into whether t z > t or t z < t and the right-hand side of Eq. (fT7|) reads, after setting x = y 

J d 3 z jf dt z [6{t - t z ) (S>G< + G<£> - S<G> - G>£<) 
+ 0& - t) (S>G> + G>£> - £>G> - G>£>)] . (21) 
The terms with t z > t all cancel, leaving 

^f^= " / dh J*dt z Tr[^>(X,z)G<(z,X)-G>(X,z)E<(z,X) 

+ G<(X,z)E>(z,X)-E<(X,z)G>(z,X)]. (22) 

This is the quantum Boltzmann equation describing the time evolution of a lepton number 
asymmetry ri£.. The initial conditions are set at t — when the lepton asymmetries and 
the RH neutrino abundances are assumed to be vanishing. All the information regarding 
lepton number violating interactions and CP- violating sources are stored in the self-energy 
S^. The kinetic Eq. (T22]) has an obvious interpretation in terms of gain and loss processes. 
What is unusual, however, is the presence of the integral over the time where the theta 
function ensures that the dynamics is causal. The equation is manifestly non-Markovian. 
Only the assumption that the relaxation timescale of the particle asymmetry is much 
longer than the timescale of the non-local kernels leads to a Markovian description. A 
further approximation, i.e. taking the upper limit of the time integral to t — > oo, leads to 
the familiar Boltzmann equation. The physical interpretation of the integral over the past 
history of the system is straightforward: it leads to the typical "memory" effects which 
are observed in quantum transport theory [22| [28| [29] . In the classical kinetic theory 
the "scattering term" does not include any integral over the past history of the system 
which is equivalent to assume that any collision in the plasma does not depend upon the 
previous ones. On the contrary, quantum distributions possess strong memory effects and 
the thermalization rate obtained from the quantum transport theory may be substantially 
longer than the one obtained from the classical kinetic theory. The very same effects play 
a fundamental role in electroweak baryogenesis 



4. The quantum Boltzmann equation for the right-handed neutrinos 

To get familiar with the out-of-equilibrium technique, we derive the quantum Boltzmann 
equation for the abundance of RH neutrinos. 

The Lagrangian we consider consists of the SM one plus three RH neutrinos N a 
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(a = 1,2,3), with Majorana masses M a . The interactions among RH neutrinos, Higgs 
doublets H, lepton doublets ^ and singlets (i = e, /i, r) are described by the Lagrangian 

= KiNJiH + hMiH c + \M a N a N a + h.c. , (23) 

with summation over repeated indeces. The Lagrangian is written in the mass eigenstate 
basis of RH neutrinos and charged leptons. 

We focus here on the dynamics of the lightest RH neutrino Ni. To find its quantum 
Boltzmann equation we start from Eq. (j7j) for the Green's function of the RH neutrino 



-M 1 )G< 1 (x,y)= - 



d 4 z 



-St (x, z)Gl (z, y) + £1 (x, z)GL (z, y) 



= j ' d 3 z J dt z [E^ 1 (x,z)G^ 1 (z,y) -E^ 1 (x,z)G^ 1 (z,y)] . 

(24) 

On the left-hand side of this equation we perform a number of operation. We first go to 
the center-of-mass coordinates and perform a Fourier transform over the spatial internal 
coordinates f. We then insert the expression in Eq. (fl4j) for the corresponding RH neutrino 
Green's function. The real part of the left-hand side of Eq. (1241) gives, after setting x = y, 
projecting onto the positive frequencies and taking the trace over the spinorial indeces 



Re 



Tr 



2 2cl>tvi 



dt 



(25) 



The self-energy of the RH neutrino is given diagrammatically in Fig. [2] (where 
indicates the generic lepton doublet in the loop) and reads 



x,y)G t ' (x,y) 



(26) 



Inserting in the right-hand side of Eq. (1241) the expressions in Eqs. fTTTl) and fTl4|) for the 
Higgs and lepton doublet Green's functions and taking the real part of it and the trace 
of the spinorial indices, we find 



dt 



dt g 



rf 3 



P 



1 



1 



(2tt) 3 2^(p) 2u H (k-p)u Nl (k) 
10 



\M{Nx — »• lH)f 



\ / 

s / 

Figure 2. One-loop self-energy of the lightest RH neutrino, 
x [f Nl {k, t){l - f e (p, t))(l + f H (k - p, t)) 

-f e (p,t)f H (k-p,t){l- f N1 {k,t))} 

x cos [(u Nl (k) - io e (p) - u H (k - p)) (t - t z )\ 

~ -2 I dt z [ ^ ^ r 1-^(^1 ^^)| 2 

7 7 (2vr)3 2^(p)2^(k-p)^(k) 1 11 ;| 
x(/, 1 (k ) t)-/f(p)0k-p)) 

x cos [(^(k) - ^(p) - w H (k - p)) (t - t z )} . (27) 

We have made the assumption that the relaxation timescale for the distribution 
functions are longer than the timescale of the non-local kernels so that they can be 
extracted out of the time integral. This allows to think the distributions as functions of the 
center-of-mass time only. We have set to zero the damping rates of the particles in Eq. ([14 ) 
and retained only those cosines giving rise to energy delta functions that can be satisfiedj. 
Under these assumptions, the distribution function may be taken out of the time integral, 
leading - at large times - to the so-called Markovian description. The kinetic equation (1271) 
has an obvious interpretation in terms of gain minus loss processes, but the retarded time 
integral and the cosine function replace the familiar energy conserving delta functions. In 
the second passage, we have also made the usual assumption that all distribution functions 
are smaller than unity and that those of the Higgs and lepton doublets are in equilibrium 
and much smaller than unity, fefn — fPfjj- Elastic scatterings are typically fast enough 
to keep kinetic equilibrium. For any distribution function we may write / = (n/n eq )/ eq , 
where n denotes the total number density. Therefore, Eq. fT27|) can be re-written as 

f For simplicity, we neglect here the fact that thermal effects may kinematically open new channels 
beyond those at zero temperature, see [2]. 
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dn Nl 
dt 



(T Nl (t))n Nl + (r JVl (t))n^ 1 , 



(r Nl (t))= I dt z I T^^r^(t), 



~ r* r d 3 k f eq ~ 

2 



(27r) 3 2Lu e 2uj H u Nl 

(2* 



where (r^^)) is the time-dependent thermal average of the Lorentz-dilated decay width. 
Integrating over large times, t — > oo, thereby replacing the cosines by energy conserving 
delta functions [29J, 

i 

dt z cos [(u Nl -UJ£- lo h ) {t - t z )} = n5 (u Nl -u t - u H ) , (29) 

we find that the two averaged rates (T^) and (r^) coincide and we recover the usual 
classical Boltzmann equation for the RH distribution function 



f d 3 k f rf 3 p 



(30) 



Taking the time interval to infinity, namely implementing Fermi's golden rule, results in 
neglecting memory effects, which in turn results only in on-shell processes contributing to 
the rate equation. The main difference between the classical and the quantum Boltzmann 
equations can be traced to memory effects and to the fact that the time evolution of the 
distribution function is non-Markovian. The memory of the past time evolution translates 
into off-shell processes. It would be certainly interesting to perform a numerical study to 
assess the impact of the memory effects onto the final baryon asymmetry. 
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5. The quantum Boltzmann equation for the lepton asymmetry and the CP 
asymmetry 

Our goal is now to compute the right-hand side of the Eq. ( |22l) describing the evolution 
of the lepton asymmetry following the CTP approach. We start with the CP asymmetry 
source term. We will see that the CP asymmetry manifests memory effects, evolves in 
time and at large times may resemble the usual CP asymmetry expression existing in the 
literature only if certain conditions are satisfied. 

As we have already mentioned, Eq. (1221) contains the information about all possible 
interesting processes for leptogenesis, e.g. AL = 1 inverse decays, AL = 2 scatterings and 
so on. To extract the CP asymmetry, we first consider the "wave" -diagram contribution 
to the lepton doublet (see Fig. [3]). From the previous discussion, we know that this 
diagram is in fact a sum of diagrams obtained assigning to the interaction points a plus 
or a minus sign in all possible manners, taking into account that, by definition, vertices 
with a minus sign must be multiplied by —1. The sum of all diagrams has to give rise 
to theta functions which ensure that the dynamics is causal. Since baryogenesis is a 
process close to equilibrium and we know that a vanishing baryon asymmetry has to be 
recovered when the RH neutrinos are in equilibrium, we can expand at linear order the 
Green's functions of the RH neutrinos by expanding their distribution functions around 
equilibrium 5f Nl = f Nl - fifc: 



where we have neglected the dependence of the distribution function on the damping 
rate. Summing all possible diagrams contributing to the loop of Fig. [3], after a long, but 
straightforward computation, we find, for instance, 



SG^ (k, t X - ty) = 5G^ (k, ty ~ t x ) 



{(tf + MO e 



4WNi (tx—ty) 



2uj Ni 




(31) 



i (£>(X, z)G<(z, X)S<(X, z)G>(z, X)) 




x 



X 



\Gf.(x 1 ,x 2 )G^(x 1 ,x 2 ) - Gf.(x 1 ,x 2 )G%(x 1 ,x 2 ) 
[Gn 2 (%2,z) - G^ 2 (x 2 ,z)] 
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Figure 3. The wave-diagram contributing to the two-loop self-energy of the lepton 
doublet. 



x [G<(z, X)G<(z, X) - G>(z, X)G>(z, X)} 

x6{z,x 2 )6{x 2 ,x 1 )6{X,z), (32) 

where we have only retained the contribution from the RH neutrino N 2 , since we will 
focus on the resonance case in which Mi and M 2 are nearly degenerate. 

Making use of the hermiticity properties in Eq. (T4J) and inserting the various 
propagators for the Higgs fields, the lepton doublets and the heavy Majorana neutrinos, 
we find that Eq. (1221) for the lepton asymmetry - as far as the CP violating source is 
concerned - is given by 



dn, 



€ 



N 



-^(t)(V Nl ) (n m -n^), 
4 

OX) 



dt 



x fdt z f Z dt 2 f 2 dhe-^-^e^ 7 -)^ [ 

JO JO Jo J {^ n > n N 1 

r d*p i - gfr) + /H q (k - P) f d\ 1 - (q) + g(k - q) 
(2tt) 3 2o;,,(p)2a;^(k - p)u Nl (k) J (2tt) 3 2%(q)2u; if (k - q)u N2 {k) 



x 



x sin (wjvi (* - *i) + (W/ + ^Jf) (*i - ^) + u N2 (t 2 - t z ) + (u>t. + Q H ) (t z - t)) 
x Tr [M x P L fM 2 f] , (33) 

where, to avoid double counting, we have not inserted the decay rates in the propagators 
of the initial and final states and, for simplicity, we have assumed that the damping rates 
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of the lepton doublets and the Higgs field are constant in time. This should be a good 
approximation as the damping rate are to be computed for momenta of order of the mass 
of the RH neutrinos. As expected from first principles, we find that the CP asymmetry 
is a function of time and its value at a given instant depends upon the previous history 
of the system. 

Performing the time integrals and retaining only those pieces which eventually give 
rise to energy- conserving delta functions in the Markovian limit (as in the previous Section, 
we do not include here the new channels that thermal effects may eventually open), we 
obtain 

3 



(t) 



dt, 



COS [(u Nl - LO ti - U H ) (t - t z )\ 



(r 2 N2 + (u N2 



X 



d 3 k g 

(2vr) 3 n e ^ 



H 



T Ni sin [(u N2 - uj Nl )t z 



)(2(u N2 - 
d 3 p 



u Nl )sm 



1-/?(P) + Z?(k-P) 



(27r)3 2^,(p)2^(k-p)^ 1 (k) 



X 



d'q l-y?(q)+/?(k-q 



Tr (M x P L fM 2 (f) . 



(34) 



(27r) 3 2%(q)2a; if (k-q)a; i v 2 (k) 
From this expression it is already manifest that the typical timescale for the building up 
of the coherent CP asymmetry depends crucially on the difference in energy of the two 
RH neutrinos. 

If we now let the upper limit of the time integral to take large values, we neglect the 
memory effects, the CP asymmetry picks contribution only from the on-shell processes. 
Taking the damping rates of the lepton doublets equal for all the flavours and the RH 
neutrinos nearly at rest with respect to the thermal bath, the CP asymmetry from the 
"wave" -diagram reads (now summing over all flavour indices) 



Im(AAt)J Mi 



x 



(AAt) n (AAt) 22 M 2 
AMt 



N 2 



2 AM sin 2 



(AM) 2 + 
rjv 2 sin [AMt] 



(35) 



where AM = (M 2 — Mi). The CP asymmetry (|35|) is resonantly enhanced when 
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H v 



Figure 4. The vertex-diagram contributing to the two-loop self-energy of the lepton 
doublet. 



AM ~ r 



tv 2 and at the resonance point it is given by 




1 Im(AAt)' 2 

2 (AAt) n (AAt) 



22 



(1 - sin [AMt] - cos [AMt]) 



(36) 



The timescale for the building up of the CP asymmetry is ~ 1/ AM. The CP asymmetry 
grows starting from a vanishing value and, for t ^> (AM)" 1 , it averages to the constant 
value quoted in the literature plj. This is true if the timescale for the other processes 
relevant for leptogenesis is larger than ~ 1/AM. In other words, one may define an 
"average" CP asymmetry 



where r p represents the typical timescale of the other processes relevant for leptogenesis, 
e.g. the AL = 1 scatterings. If r p ^> 1/AM ~ T^, the oscillating functions in ( |36i) 
are averaged to zero and the average CP asymmetry is given by the value used in the 
literature. However, the expression ( 135]) should be used when^ p l/AM ~ . 

The fact that the CP asymmetry is a function of time is particularly relevant in 
the case in which the asymmetry is generated by the decays of two heavy states which 
are nearly degenerate in mass and oscillate into one another with a timescale given by 
the inverse of the mass difference. This is the case of resonant leptogenesis [31] and soft 
leptogenesis [32J . From Eq. fl35l) it is manifest that the CP asymmetry itself oscillates with 
the very same timescale and such a dependence may or may not be neglected depending 
upon the rates of the other processes in the plasma. If > Tn 2 , the time dependence 
of the CP asymmetry may not be neglected. 




(37) 
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Figure 5. One-loop diagram contributing to the self-energy of the lepton doublet. 



The expression (j3"5|) can also be used, once it is divided by a factor 2 (because in the 
wave diagram also the charged states of Higgs and lepton doublets may propagate) and 
the limit M 2 3> Mi is taken, for the CP asymmetry contribution from the vertex diagram 
(see Fig. BJ 



(*) 



Im (AAt)J Mi 



2 sin 2 



Mot 



1 N 2 

M 2 



sin [M 2 t] 



(3? 



167r(AAt) u M 2 

The timescale for this CP asymmetry is ~ M2 and much larger than any other timescale 
in the dynamics. Therefore, one can safely average over many oscillations, getting the 
expression present in the literature 



(el) 



Im(AAt)' Mi 



(39) 



167r(AAt) n M 2 ' 

Finally, the AL = 1 inverse decays can be computed along the same lines described 
previously by considering the one-loop contribution to the lepton doublet self-energy E| 4 
(see Fig. The equation for the lepton asymmetry then becomes 



eq 

(n Nl -nZ) ~ (T™(t))^n Ci (40) 



^ = 4(*)<r*> K - r#J - <r£(t)> 



where eV (t) = e]y* + eff^ is the total time-dependent CP asymmetry for the flavour i, and 



(T Ni (t))-2l dt,, (27r)3 ^ y (27r)3 2^2^^ 

x cos^-w^-Wff)^-^], (41) 
is the time-dependent thermal average of the inverse-decay interaction rate. 
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6. Conclusions 

The quantum Boltzmann equations derived in this paper can be used to perform a 
thorough investigation of the impact of flavour effects onto leptogenesis and, on more 
general grounds, to provide a quantitative relation between the light neutrino properties 
and the final baryon asymmetry. It would be interesting to see how large are the 
corrections to the baryon asymmetry once the non-Markovian description is adopted, 
including memory effects and off-shell corrections. They may lead to the slowdown of 
the relaxation processes thus keeping the system out of equilibrium for longer times and 
therefore to an enhancement of the final baryon asymmetry. It would also be of interest 
to see the impact of our results on the transition between one flavour and two flavours as 
we discussed in the Introduction. 

One of the main results of our investigation is that the CP asymmetry turns out to 
be a function of time and its value at a given instant of time depends on the past history 
of the system, see for instance Eq. fl34l) . This result is relevant when the timescale of the 
evolution of the CP asymmetry is larger than the timescale of the other processes. We 
have pointed out that this is relevant when the asymmetry is generated by the decays of 
two nearly mass-degenerate heavy states and the resonant effects are exploited. 
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